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Abstract 

A method for the calculation of the conductance of nanoscale electrical junctions is extended 
to ab-initio electronic structure methods which make use of the periodic supercell technique, and 
applied to realistic models of metallic wires and break-junctions of sodium and gold. The method 
is systematically controllable and convergeable, and can be straightforwardly extended to include 
more complex processes and interactions. Important issues about the order in which are taken 
both the thermodynamic and the static (small field) limits are clarified, and characterized further 
through comparisons to model systems. 

PACS numbers: 
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I. INTRODUCTION 



Nanoscale and molecular electronics is one of the most active topics of research in physics 
today™ 1 ; it has very important consequences, both for fundamental research and for in- 
dustrial processes, which will reach their inherent quantum limits in the decade to come. 
Accurate methods of theoretical as well as experimental characterization are essential, and 
many open questions remain about the structure, equilibrium, and dynamics of nanometer- 
sized systems carrying electrical currents. The ab initiosimulation of materials properties 
is in a unique position to develop the scope and our understanding of electronics at the 
nanoscale, providing the only method of systematic analysis of electronic and structural 
effects, which are never all simultaneously accessible in experiment. Present simulations of 
nanoscale transport usually fall into two categories, either employing time-dependent den- 
sity functional theory (TDDFT)3EE1 

or various flavors of Landauer-Biittiker-like formulas^^, 
sometimes using a non-equilibrium Green's function formalism^! (for a review see Ref. |9]). 
Due to the way they are formulated, using embedding schemes, these techniques often rely 
on localized functions (atomk P 11 * 12 * 13 * 14 -! or Wannier^HSl functions) to describe the single 
particle wavefunctions of the system. This introduces an inherent difficulty in converging 
calculations, as the basis sets usually cannot be refined systematically. There is still a great 
deal of uncertainty about the precision of both experiments, e.g. due to fluctuation in exper- 
imental conditions of contact and current flow, and theory, where no standard model is yet 
accepted as being predictive of experiments, apart from simple cases of continuous contact 
with conductances of at least one quantum of conductance (Go)- The role and importance of 
electron-electron 1 9 1 17 1 18 1 19 1 2 -*^ 1 and electron- vibration 1 22 1 23 1 24 1 interactions, has been recognized 
as an important factor in obtaining the correct transport properties even though satisfactory 
treatment for a general system at the ab initio level is not yet available. 

In the following, we present a method to calculate the conductance of a quantum junction, 
which can be systematically converged and extended to include the effects of different inter- 
actions. The formalism has been previously applied to model systems in lD 2 ^ and jellium 
slabs 2 -^, and is here extended to incorporate 3D realistic ab initio electronic structures at the 
level of local or semi-local time-dependent density-functional theory. Section [TT] describes the 
method and how it must be adapted to suit the periodic boundary conditions and supercells 
which are often used in ab initio calculations. Preliminary numerical characterization is 
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carried out in Section III Section IV analyzes the convergence properties of the method and 
applies it to monatomic wires of sodium. In particular, the possibility of obtaining precise 
calculations of very low conductances is explored for a tunnel junction. Finally, Section [V] 
examines contact geometry and bonding effects for bulk gold electrodes. 

II. METHODOLOGY 

The conductance of a nanoj unction characterizes the long-time dynamics of the electronic 
response of electrons to a driving electric fielcP^. While formally the long-time limit demands 
the study of an extended system, for calculation purposes it is possible to consider a finite 
model for a finite time and the resulting conductance is obtained by extrapolation of the 
conductance function 

G 2P = lim lim G 2P (oo,L) } (1) 

uj^0+ L—>oo 

where G2P is the two-point conductance (see RefP*j). The order of limits is important here: 
the one given characterizes transport in an extended system where as the reverse would 
reflect damped oscillations of density in a finite (even though large) system. The underlying 
finite system can fulfill any boundary conditions and these do not affect the extrapolated 
results. We take advantage of this fact and use the periodic boundary conditions and a 
plane-wave basis for the ab initio calculations below. 

It is presumed in the following that the junction whose conductance we are searching for 
is centered at in the middle of the cell (which thus extends from —L/2 to +L/2). 

The calculation of the electronic response function at the level of local or semi-local TD 
DFT, which then leads to the conduction function (Eq. [I]), proceeds in three steps. First 
we perform a calculation of the occupied (e n < 0) as well as the unoccupied (e n > 0) 
eigenenergies e n and eigenstates n (r) of the system. 

Second, we use the eigenvalues and eigenstates to compose the positive imaginary-time 
Matsubara Green's function 

G(ry;r) = J2MrW n (r')^ r ^- 1 (2) 

n 

As Q is anti-periodic in imaginary time (fermionic) there is no need to specify explicitly its 
behavior for negative imaginary time. 
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The electronic response to the total electric field is characterized by the polarizability 

P(r,r';r) = -Q(r,r';r)G{r' ,r;-r) 

= g(r,r';r)g(r',r;(3-r) (3) 

which, after Fourier transformation, r — > u, and integration over the cross-sectional area of 
the junction, A, gives the integrated polarizability relevant for charge transport 

P(x,x';iu) = / / ' dS±dS' ± P(r,r';iu). (4) 



Methods going beyond the present level of approximations, i.e. using nonlocal exchange- 
correlation kernels^^ or Green's-function-based many-body methods^, would differ in the 
above Equations [2] and [3] The expression for the irreducible polarizability, Eq. [3j would 
contain further vertex diagrams^ and, in the case of many-body methods, the Green's 
function cannot be expressed in terms of one-electron wavefunctions as in Eq. [2} However, 
the discussion that follows would apply also to these computationally more demanding 
approaches. 

Finally, the third step consists of integrating the polarizability to obtain the conductance 
function. For an infinitely long system, the conductance is obtained from the expression^! 

i2P i 



G (u>, oo) = uj I I P(x,x';iu>)dxdx', (5) 

J -oo JO 

where this integral converges for any finite u since P(x,x'; iou) — > for x, x' — > oo. The in- 
tegration region corresponds to choosing elements of the polarizability which connect points 
on opposite sides of the junction. This is intuitive, as we are interested in how a perturbation 
on one side can influence charges on the other side, through the junction. 
For a finite system of length L we obtain the corresponding function 

G 2P (lu,L)=lu / P(x,x';iu)dxdx', (6) 
Jv 

where T> is a domain of positive x' and negative x which must guarantee the correct limiting 
procedure. Using a periodic supercell, for x 1 — x — > ±L, we approach a periodic image 
of the system we want to study. Further, if the system is not translationally invariant we 
will approach a region of the (x, x') plane (the lower right hand corner of Fig. [T]) where 
the polarizability behaves very differently from that near (0, 0) (typically one with a more 
metallic behavior and larger polarizability than the junction). Hence, a correct integration 
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FIG. 1: Left: example of a color plot of the polarizability P(x, x') for a junction. P is non zero only 
near the diagonal x = x' , and lower in the central tunneling region (see below). Right: the region of 
spatial integration for the polarizability as in Eq.[6} For periodic systems the polarizability will have 
spurious images, which must be excluded from the integration region. Taking the thermodynamic 
limit will increase the size L of the system, and the triangular domain converges uniformly to the 
quarter plane < x < oo and — oo < x' < 0. The inset is a cartoon of an atomic wire with a 
central site and two positions, x and x'. 

needs to truncate the quadrant denned by —L/2 < x < and < x' < L/2. For a finite 
system there is no unique choice of T>, but there is a natural one, which is < x' < L/2 
and —\x'\ < x < 0, or equivalently < x' — x < L/2, defining a triangle between and the 
points (0,L/2) and {-L/2, 0) (Fig. [I}. 

The finite size of the system determines the minimal frequency which can be reliably 
described in the conductance, or equivalently the longest time propagation. For longer times 
or lower frequencies the electrons will reach the limits of the system, and the conductance 
decays. The minimum frequency may be estimated as: 

u min = 2-Kv F /L, (7) 

where vf is the Fermi speed. Thus, an electron at the Fermi level takes time l/u m i n to 
traverse the whole system. This frequency will be essential in determining how to extrapolate 
the conductance function to zero frequency. 

We apply the formalism described above to systems described by modern electronic struc- 
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ture methods. Many of these techniques use periodic boundary conditions to describe crys- 
talline structures, and represent wavefunctions and electronic densities using plane wave 
basis sets. Here we describe the corresponding small changes needed in the formalism. 

First, as we wish to describe an isolated nanojunction between two leads (which are in 
principle infinite) we will use only the zone-center T k-point of the Brillouin zone (BZ) along 
the axis of conduction. Using several k-points would in effect simulate an array of interfering 
junctions. The periodic boundary conditions are nevertheless exploited as the regions near 
the edge of the simulation cell are described continuously instead of being brutally cut off or 
terminated with hydrogen atoms. The thermodynamic limit along the conduction axis must 
still be ensured by increasing the longitudinal system size until the conductance converges. 

In the directions perpendicular to x a denser k-point grid can be used if bulk three- 
dimensional leads are considered: the electronic structure of the leads will thus be repre- 
sented correctly, but care must be taken that the transverse distance between images of the 
"junction" part of the cell is sufficient to avoid interference between periodic images. In 
the case of a purely ID system no perpendicular k-points are necessary as the system is 
supposed to be isolated in vacuum along y and z. 

III. NUMERICAL CHARACTERIZATION 

In order to understand the convergence behavior of the main results we have also analyzed 
a finite ID jellium model and a finite ID tight-binding model. The length L of both systems 
can be made much larger than in the ab initio models, which allows detailed study of the 
extrapolation to small frequencies. Similarly to the ab initio case, both models use periodic 
boundary conditions and their parameters are such that the density and the Fermi speed 
of the particles will be identical to that of a sodium chain studied within the self-consistent 
ab initio calculations. In fact, the jellium model with a bare electron mass is an excellent 
model for the sodium wires. This is shown in the Fig. |2]where the dispersion of eigenenergies, 
obtained from the ab initio calculations, is compared with the dispersion of the ID jellium 
and the fitted TB model. More generally, these two models represent two extremal types of 
electronic structure for metallic wires. The correct understanding of the extrapolations to 
infinite size and zero frequency of their conductance functions is very useful for performing 
extrapolations of more realistic but numerically more demanding ab initio calculations. 
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FIG. 2: The dispersion of eigenenergies, obtained from the ab initio calculations (see Sec. IV for 



details), is almost identical to that of ID jellium using the bare electron mass. The TB model is 
fit to have the same Fermi speed, and the differences with respect to the ab initio dispersion are 
more significant. 

For both model systems one can find the eigenstates exactly by going into reciprocal 
space. The conductance function (Eq. [6]) can be expressed using the exact eigenstates of the 
Hamiltonian 

G 2P (uj,L) = 2u^2 \ s (Pi - Qj)\ 2 { 1 -n qj )n Pi 

e P{e p .-e qj ) _ j 



i.l 



X 



iuo + e 



(8) 

where the sum goes over all eigenstates, Pi or ^ are the momenta of the eigenstates, n n 
are Fermi occupancies of the state Pi at temperature T = 1/ksP, e Pi are the eigenenergies 
corresponding to an eigenstate with the momentum pi, and the factor 2 accounts for the 
spin degeneracy. The function s( ) represents the conductance vertex- factor (similarly to 
the expression for the conductivity in terms of the polarization function, see e.g. Bruus and 
FlensbergSS) and takes different forms for the two systems, as given below. 

The jellium model consists of a ID non-interacting electron gas of total length L and 
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density n = N/L where N is the total number of electrons. The eigenvalues are 

Pi 27T , , 

e K = ^-E F , Pi =— m = ±1,±2,... (9) 

and the conductance vertex-factor is 

s(p) = — - / dxe~ ipx = p=— (1 - e - ipL/2 ) . (10) 

The Fermi energy is obtained from the requirement that the charge per unit length is iden- 
tical to that of the sodium wires; the length L is set to Nd^ a , where d^ a is the inter-atomic 
distance of the sodium wire so that the density of electrons is n = N/L = l/d^a- 
In the case of the tight-binding model, the Hamiltonian has the form 

N/2-1 

H TB = £ -'-(clc^ + h.c), (11) 

-N/2 

with the resulting eigenvalues e n = —tcos(pi) for a state with momentum pj. The momen- 
tum only takes discrete values 

2tt 

Pl = — i, i = 0,l,2,...,JV-l (12) 



and the conductance vertex-factor is 

I ^11 | MN/2+1) _ i 

■<P) = e ipn = A=~ = • (13) 



V2N ^ v/2 - 1 



At half filling, the Fermi momentum fcp = P(n/4) and the Fermi speed is t>i? = de(p)/dp = t. 



Making use of the Fermi speed of the sodium wires considered in Sec. [IV]t>F = 0.33 a.u., we 
identify the TB parameter as t = VF/d^a ~ 0.07. 

The evaluation of the expression ^ for both models is very fast and can be done for 
much longer wires than in the case of first-principles calculations. In the top panel of Fig. |3j 
we show the conductance functions for jellium wires of lengths L = Nd^ a with N = 4, 8, 16 
and 100 at temperatures much lower than the Fermi energy. The curves converge smoothly 
to the zero-temperature infinite-length limit that is known analytically^, and readily give 
the static limit of one quantum of conductance G 2P (0,0) = 2e 2 /h = 1/tt a.u. Furthermore, 
the functional form shows finite size effects precisely according to the expected criterion 
(Eq.|7|) 

44 

Umin = = 0.11, 0.05, 0.03, 0.004 (14) 
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FIG. 3: Convergence of the conductance of a jellium wire with respect to system size. The lengths 
correspond to sodium wires 4, 8, 16 and 100 atoms long. Upper graph: T <C Ep {(3 = 1000). Lower 
graph: T ~ Ep {(3 = 30). The curves at low temperature approach the analytical infinite-size result 
(for T=0 - continuous line). The curves should be extrapolated to zero frequency, disregarding 
values for u < uj m in (indicated by a dot on each curve) . 
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for iV = 4, 8, 16, 100 respectively. The extrapolation for N = 8 or N = 16 gives good 
estimates of the conductance, through the limit in Eq. (|6|. At temperatures comparable 
with the Fermi energy (bottom graph in Fig. [3]) the extrapolated value is somewhat below 
the zero temperature limit but the functional form of the conductance is essentially identical. 

On the other hand, the tight-binding model, shown in the Fig. [4], offers less reliable ex- 
trapolations. This is caused by the non-monotonic behavior of the conductance function 
at small frequencies, which in turn arises because of the bandwidth of the model (here the 
bandwidth is 2t = 0.14a.u.). This non-monotonic behavior weakens if we look at the conduc- 
tance curve at higher electronic temperature. This is shown in the lower graph in the Fig Q 
for a temperature comparable to the bandwidth. For these temperatures the conductance 
function can be easily extrapolated to (high-T) conductance values identical to the jellium 
model (Fig. |3j lower graph). This suggests that, in principle, by performing the calcula- 
tions at different electronic temperatures or smearing, one may enhance the extrapolation 
procedure in realistic calculations, which may combine aspects of free-electron (jellium) and 
localized electron (TB) behaviors. 

To summarize, the models show good convergence in the extrapolated conductance at 
uj = for systems of length equivalent to 8 or 16 atoms. The jellium dispersion is quite 
close to the ab initio one for Na, whereas the tight binding one is not - this could be 
expected from the simple metal nature of sodium. The conductance function of the TB 
chain is qualitatively different, and actually overshoots the quantum of conductance for 
small imaginary frequencies and large L. Finally, for high temperature (T comparable to 
Ep) both models depart from the analytic curves for T=0 and the conductance decreases. 

IV. SODIUM MONOWIRES: SIZE CONVERGENCE AND ENERGY DEPEN- 
DENCIES 

We begin the ab initio studies with a prototypical application: the calculation of the 
conductance of a uniform monatomic wire (monowire) of sodium atoms. With one s electron 
per atom, and given the simple-metal nature of sodium, the conductance in the independent 
particle case will be two quanta of conductance, due to spin-degeneracy. 

Technical details 

The ground state wavefunctions and electronic structure are calculated within density 
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FIG. 4: Convergence of the conductance function of a tight-binding model with respect to system 
size (4, 8, 16 and 100 atoms), vf is fixed to that of the sodium wire. Upper graph: T <C It 
{(3 = 1000) Lower graph: T ~ 2t (/5 = 30) where 2t is the bandwidth. The curves give the correct 
zero temperature conductance for L — > oo, but the extrapolation to u = is complicated by the 
non-monotonic behavior at small frequencies. 
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FIG. 5: Convergence of the conductance of a continuous monowire of sodium atoms, with respect 
to cell size (4, 8, and 16 atoms in the unit cell). The minimum frequency for which the conductance 
function is valid depends on the Fermi speed and goes down with increasing system size. For 8 
atoms the extrapolation to zero frequency is already good, arriving close to the expected 2 quanta 
of conductance (horizontal line) for a non-interacting system with spin-degeneracy. 



functional theory (DFTp 11221 , using a plane wave representation, with the ABINIT^ or 
SFHINGXEl codes (the results have been checked to be independent of the ground state 
code used). We employ norm-conserving Troullier-Martins^ type pseudopotentials, with 
nonlinear core corrections^ and the d channel as a local potential. The kinetic energy cutoff 
(20 Hartree) and number of bands (100 per Na atom) were over-converged to allow for full 
checks of the convergence of the conductance calculation. The calculation of the conduc- 
tance, in a module of the GWST code^, was carried out with a kinetic energy cutoff of 8 
Ha. 

The inter-atomic distance was set to 2.477 A. Other distances were checked, but do not 
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influence the results appreciably: the Fermi point for the wire is fixed by the parabolic nature 
of the bands, and, most importantly in our case, the Fermi speed scales as the inter-atomic 
spacing, making the critical minimum frequency (Eq. [7]) independent of the spacing. 

The perpendicular size of the unit cell is more important, as we wish to simulate a truly 
ID system using a supercell. The lateral dimensions of the unit cell are fixed to 4 A, which 
is enough to ensure that conductance only happens along the wire direction. Checks with 8 
and 16 A cells showed that the conductance function is already well reproduced to within a 
few percent. 

Results for uniform Na monowires 
Figure [5] shows the size convergence of the conductance function, for unit cells containing 
4, 8, and 16 atoms. As can be seen in this simple case the conductance function is correct 
to lower and lower frequencies as the system size is increased, and the extrapolation tends 
towards 1/ir. Already with an 8 atom cell the linear extrapolation of G 2P to zero frequency 
gives a value very close to l/V, as expected from the previous section. Comparing the 
values obtained for successively larger system sizes gives an estimate of the residual error 
in converging L. A 16 atom unit cell is about 40 A long. The very slow convergence of 
G(u m i n , Ij) as a function of L is due to the ID character of the system. Coulomb screening 
in 1 dimension is quite inefficient and the polarizability decays quite slowly. In 3D systems 
the screening will be stronger and the size convergence quicker 

A very important point is to have a good estimate of the Fermi level. In ID atomic 
chains this is not trivial, as the equivalent of k-point sampling is the length of the system. 
The distance between the levels bracketing the Ep decreases with L, but not uniformly. For 
small L, there is an alternation in the position of the Fermi level: for wire lengths which are 
multiples of 4 the Fermi level is exactly on a single particle state, whereas for other values it 
is higher, and between states. The dielectric response of these two cases is very different, as 
one case appears to be a metal and the other an insulator (whose gap goes to as L —>■ oo). 
In the interest of brevity we have not included the results for L = Ndpf a with N = 5, 6, 7 in 
Fig. [5j they oscillate slightly (as a function of N) and converge more slowly, though to the 
same end result. We will see these effects again in the next Section for the case of a wire 
with a gap. 

In our formulation the conductance is expressed in imaginary frequency. This implies 
contributions from all electron hole pairs in the polarizability, not just those for states near 



13 




FIG. 6: Convergence of the conductance of a wire with a gap (1 atom missing), with respect to 
the size of the lead wires (2 to 8 atoms in each lead). The conductance functions converge well for 
even numbers of atoms in the leads, due to a correct positioning of the Fermi level, and can be 
extrapolated to zero frequency (the line segment extrapolates from 8-atom lead case). With odd 
numbers of atoms the Fermi level is positioned in an artificial gap. 

the Ep. Because of this the convergence in the number of states is comparable to (but 
slightly faster than) that of a GW calculation, with between 5 and 20 bands per electron 
(for Na wires we need 10 bands per electron). 

Na monowires with a gap 
We now proceed with an inhomogeneous case: a wire of Na atoms with a gap (of width d). 
This is the simplest example of a nanojunction. The conductance will naturally go into a 
tunneling regime as the gap becomes wider. This example is important firstly because it 
has a simple dielectric response, but also because tunneling is an important and extreme 
regime for the conductance. Fig. [6] shows the size convergence of a wire with a gap of one 
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FIG. 7: Conductance function of gapped wires, with respect to imaginary frequency, and for 
different gap sizes (in fractions of the inter-atomic distance in the regular wire). There are 8 
atoms in each lead. Inset: extrapolated frequency conductance, as a function of gap length. The 
conductance saturates for small gaps, and for gaps larger than ~ 0.75 inter-atomic distances the 
tunneling decay appears. 

atom. The two remaining parts of the wire are of equal length, increasing from 2 to 8 
atoms (each). The even length leads converge relatively quickly to a regular conductance 
function, from below. The odd length leads converge from above but very slowly, because 
the HOMO/LUMO states are quite far apart, which gives a badly placed Fermi level, as 
above for continuous wires (bad in the sense that it is far from the limiting Fermi level for 
an infinite system). Extrapolation to frequency gives a conductance of 0.05 a.u. (±0.005). 

The oscillations of G{uS) are due to aliasing effects in the Fourier transform from imaginary 
time to imaginary frequency. These are in some cases difficult to eliminate for very low 
amplitude elements of P(x,x' ,ir). 
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Finally we consider the transition from the conducting to the tunneling regime, by in- 
creasing the gap size (for 8 atom leads). The conductance functions are represented in Fig. [7j 
with the inset showing the decrease of the extrapolated static conductance as a function of 
the size of the gap. The conductance decreases with gap size, and by a gap of 3/4 of an 
atomic spacing (about 4 Bohr) the tunneling regime begins with exponential decay of G. 
The tail of G(d) gives a very good fit to G(d) = 1.68 exp(— 0.770 d) where d and G are in 
atomic units. The dashed lines are linear extrapolations of G fit to the interval [0.05, 0.2] 
Hartree. For the lowest curves (largest gaps) the fit was performed further out, on [0.2, 1], 
as the functions are flatter and the aliasing noise more important. In this way we are able 
to represent quite small conductance values, down to 0.001 a.u. (or 0.003 Go). 

A similar system was examined by Beste et al. in Reference HH but with gold chains 
instead of Na. They find a conductance of about 0.12 a.u. for a gap of d^u/2 (1.28A), 
which is close to our value of 0.18 a.u. for d^ a j2. One important conclusion of Ref.^is the 
very strong deviations which can appear depending on the basis set nature, with localized 
basis sets. Our results can be converged systematically, using a plane- wave basis set, but 
are probably heavier calculations as a trade-off. 



V. GOLD WIRES: LEAD STRUCTURE AND K-POINT SAMPLING 

We now proceed with a more structured system, showing an explicit constriction. A gold 
junction is made from a 2 atom wire contacted to bulk 3D electrodes. The electrodes are 
FCC stacked gold (at the experimental nearest neig hbor distance of 2.9 A of Ref. 138} of 
which we use a 2x2 (111) surface unit cell. The wire atoms are evenly spaced with the FCC 
inter-layer distance of 2.37 A (which is compressed compared to the DFT-LDA equilibrium 
distance of 2.5 A for the infinite straight wire®'). No relaxation of atomic positions with 
respect to the bulk is taken into account, but the addition of further structural effects is 
in no way more difficult; contrary to some other approaches to transport 4 14 we are not 
constrained to specific unit cell lengths or layer spacings. A typical unit cell is shown in 
Fig. [8j for the minimal electrode thickness of 2 layers (in each electrode). Electrodes 3 and 4 
layers thick were also tested. Because of periodicity, and in order to maintain a continuous 
FCC structure at the cell boundary, the point of contact of the wire to the right electrode 
alternates between the different possible FCC stacking sites (for 2 and 4 layer electrodes), 
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FIG. 8: Unit cell of a short gold wire contacted with bulk gold electrodes, which are 2 layers thick. 
Both atoms sit at natural FCC hollow positions on the (111) surfaces. 

and an on-top position (for the 3 layer electrodes). We find little effect of contact position 
on the conductance of the junction (see below) in this continuously metallic, well contacted 
case. 

Technical details 

The pseudopotential we use is of the Hartwigsen-Goedecker-Hutter flavor^!, with only the 
6s electrons in the valence. Our choice of pseudopotential is justified by its softness, by 
the chemical homogeneity of the system, and our intention to go beyond LDA and include 
many-body corrections. We have performed tests on FCC gold in the GW approximation 
(which is beyond the scope of the present article), including the 5d electrons. The problems 
recognized by Marini et al^ for Cu appear for Au as well: the exchange self-energy is quite 
badly represented for the 5d electron states, due to the absence of the 5s and 5p. The latter 
are far in energy but have an important spatial overlap with the 5d. Consequently, the 
exchange self-energy lacks important contributions if one uses only the valence electrons. 
The GW d bands are very poor (whereas their position in LDA is very close to experimental 
values); some bands are pushed down and others up to the Fermi level, which would change 
the conductance severely. The use of a purely 6s potential is less realistic but reduces these 
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exchange effects (which are now between the 6s and the core 5d states). A more complete 
solution is that adopted by Shishkin and Kresse^ in the PAW formalism 43 . As PAW allows 
explicit reconstruction of the core states, the exchange with the valence can be calculated 
explicitly. Finally, the d electrons do not complexify the independent particle conductance 
calculation formally, but do make the calculations much heavier (with an additional 10 elec- 
trons per atom). As the states at the Fermi level are purely s-electron like, the conductance 
will not be affected strongly. However, as our method is an integral of the dielectric response 
of the system, the absence of the d electrons will have an indirect effect, through changes in 
the polarizability. 
Results 

From a calculation of a uniform wire (with k-points along the wire axis), we estimate the 
Fermi speed in the wire to be 0.42 au (1.9 10 6 m/s), corresponding to a wavelength of 7.22 
Bohr. A simple metal approximation for the bulk gives an estimated Fermi speed of 0.64 
au (1.4 10 6 m/s) from the Au Seitz radius. A DFT calculation of FCC bulk naturally gives 
a more complex band structure - the modulus of the Fermi speed varies by some 20% in 
reciprocal space. The HGH pseudopotential gives an average value of 1.02 au (2.2 10 6 m/s). 
A more complete pseudopotential with d electrons reduces this value to 0.67 au (1.5 10 6 
m/s). The value we are interested in is the speed of propagation of an electronic signal 
through the whole system, i.e. through the 3D bulk (with the pseudopotential we are using) 
and the wire, which will be between the pure bulk and pure wire values. With the bulk and 
wire Vp, we can estimate the minimal frequencies which can be represented for different unit 
cell sizes. With cells of lengths 27, 36, and 45 Bohr, we obtain u min = 0.163, 0.132, and 
0.113 Hartree. 

In Fig. [9] (lowest three curves) we represent the conductance as a function of imaginary 
frequency for a series of gold junctions (like that schematized in Fig. [8]) with 2, 3 and 4 
layers of FCC gold in the bulk leads. As before for linear wires, initially only the T point 
wavefunctions are used, and G(uj) takes similar values to the case without leads, which would 
suggest a similar system length of ca. 8 layers in each lead to converge the conductance. The 
extrapolated value, close to the quantum of conductance is in agreement with the results of 
the more extended models of the monatomic gold contacts^. 

With bulk leads it is essential to look at k-point convergence: to represent the bulk states 
correctly, we increase the sampling of k-points in the direction perpendicular to the junction 
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FIG. 9: Conductance functions of gold junctions containing a 2 atom wire and bulk FCC leads. 
For 2, 3 and 4 layers of gold in the leads, and only the T k-point (solid black, dotted, and dashed), 
and for 2 layers of gold and a 4x4 sampling of the BZ perpendicular to the wire axis (solid orange 
curve). The horizontal line is the quantum of conductance. The peak in the conductance moves 
to lower energies as the system length is increased. In the case with denser BZ sampling the 
conductance function is much better represented, even if the minimum frequency is not changed. 

axis. Again, as noted above, one must also take care to keep the nanojunctions themselves 
well isolated in the perpendicular direction, to avoid interference effects which would be am- 
plified by the perpendicular k-points. With 2 layer leads and a 4x4 sampling of k-points, we 
find the conductance function shown by the top curve of Fig. [9j The low frequency behavior 
is now much closer to that of the longer Na wires. Increasing the lead size thus has two 
effects, which are controllable separately in this 3D case: first improving the representation 
of the density of states (DOS) (which can also be achieved by using the perpendicular k- 
points if the leads are bulk-like), and second lowering the minimum representable frequency 
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Umin (which can only be achieved by increasing the system length L). In the future, a more 
extended study will combine the effects of longer leads and perpendicular k-point sampling, 
but the computational load will require the parallelization of our code, which is underway. 

To summarize, the bulk leads on gold wires show that even with relatively small system 
sizes a constriction limits the conductance to a single quantum of conductance. Bulk 3D 
leads give much stronger screening than ID ones, and faster convergence of the conduc- 
tance function. Including k-points to sample the perpendicular electronic states in the leads 
improves the description of the DOS and the screening. 

VI. CONCLUSIONS 

We present a new method to calculate the transport properties of nanoscopic junctions. 
The extension of the basic formalism to periodic systems and realistic electronic structure 
is detailed, and the convergence properties are compared to model systems. Order-of-limits 
problems are reviewed, which also concern many other approaches to quantum transport, 
as well as numerical issues. The differences and inherent advantages of the method are 
discussed, in particular the way it treats leads and its systematic convergeability in number 
of single particle states and the spatial representation of different quantities. Applications 
to sodium and gold nanojunctions is presented. The first show the properties of purely ID 
systems, and demonstrate the variation between regimes of continuous metallicity and of 
tunnel junctions. The gold junctions explore contact geometries and some of the fundamental 
differences between ID and 3D electrode structure and screening. 
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